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Abstract 

We present exact solutions to the Einstein-Maxwell system of equations with a 
specified form of the electric field intensity by assuming that the hypersurface 
{t = constant} are spheroidal. The solution of the Einstein-Maxwell system is 
reduced to a recurrence relation with variable rational coefficients which can be 
solved in general using mathematical induction. New classes of solutions of lin- 
early independent functions are obtained by restricting the spheroidal parameter 
K and the electric field intensity parameter a. Consequently it is possible to 
find exact solutions in terms of elementary functions, namely polynomials and 
algebraic functions. Our result contains models found previously including the 
superdense Tikekar neutron star model [R. Tikekar, J. Math. Phys. 31, 2454 
(1990)] when K = —7 and a = 0. Our class of charged spheroidal models gener- 
alise the uncharged isotropic Maharaj and Leach solutions [S. D. Maharaj and P. 
G. L. Leach, J. Math. Phys. 37, 430 (1996)]. In particular, we find an explicit 
relationship directly relating the spheroidal parameter K to the electromagnetic 
field. 



1 Introduction 

Exact solutions of the Einstein-Maxwell field equations are of crucial importance in 
relativistic astrophysics. These solutions may be utilised to model a charged relativis- 
tic star as they are matchable to the Reissner-Nordstrom exterior at the boundary. A 
recent review of Einstein- Maxwell solutions is given by Ivanov^. It is interesting to 
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observe that, in the presence of charge, the gravitational collapse of a spherically sym- 
metric distribution of matter to a point singularity may be avoided. In this situation 
the gravitational attraction is counterbalanced by the repulsive Coulombian force in 
addition to the pressure gradient. Einstein-Maxwell solutions are also important in 
studies involving the cosmic censorship hypothesis and the formation of naked singu- 
larities. The presence of charge affects values for redshifts, luminosities and maximum 
mass for stars. Consequently the Einstein-Maxwell system, for a charged star, has 
attracted considerable attention in various investigations. 

In an attempt to generate exact solutions in some cases Vaidya and Tikekar- pro- 
posed that the geometry of the spacelike hypersurfaces generated by {t = constant} 
are that of the 3-spheroid. This spheroidal condition provides a clear geometrical in- 
terpretation which is not the case in many other exact solutions. Knutsen- was the 
first to consider the pressure gradients of stars with spheroidal geometry and showed 
that they were negative. Also note that spheroidal geometries exhibit the important 
physical feature of being stable with respect to radial pulsations^. Tikekar- compre- 
hensively studied a particular spheroidal geometry and showed that it could be applied 
to superdence neutron stars with densities in the range of 10^^ g cm~^. Maharaj and 
Leach- found all spheroidal solutions, for uncharged stars, that could be expressed in 
terms of elementary functions. Mukherjee et al— showed that it was possible to express 
the general solution in terms of Gegenbauer functions, an alternate form of the general 
solution was found by Gupta and Jasim-. These uncharged solutions can be extended 
to models in the presence of electromagnetic field. Spheroidal models in the presence of 
an electric field have been extensively studied by Sharma et al&, Patel and Koppar— , 
Patel et al^^, Tikekar and Singh—, and Gupta and Kumar—. These investigations have 
been motivated on the grounds that by restricting the geometry of the hypersurfaces {t 
= constant} to be spheroidal produces neutral and charged stars which are consistent 
with observations for dense astronomical objects. Models with spheroidal geometry can 
be directly related to particular physical situations: the maximum mass is in agreement 
with values for cold compact objects^, values for densities are consistent with strange 
matter—, the equation of state is consistent with a compact X-ray binary pulsar Her 
X-1— , relevance to equation of state for stars compared of quark-diquark mixtures in 
equilibrium^^, and uniform charged dust in equilibrium—. Spheroidals geometries are 
relevant in core-envelope steller models, core consisting of isotropic fiuid and envelope 
with anisotropic fiuid, as shown by Thomas et al^^ and Tikekar and Thomas^. 

It is clear that stars with spheroidal geometry have a number of different physical 
applications and therefore require deeper investigation. In this paper our objective 
is to generate a class of charged spheroidal solutions corresponding to a physically 
reasonable form for the electric field intensity. Our intention is to obtain simple forms 
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for the solution that highhghts the role of the spheroidal parameter. In section 2 we 
express the Einstein-Maxwell field equations for the static spherically symmetric line 
element. The condition of pressure isotropy becomes a second order linear differential 
equation which facilitates the integration process. We assume a solution in a series 
form which yields recurrence relations, which we manage to solve from first principles in 
section 3. It is then possible to exhibit exact solutions to the Einstein-Maxwell system. 
In section 4 we present polynomials as first solutions and product of polynomials and 
algebraic functions as second linearly independent solutions. In addition we express 
the general solutions in terms of elementary functions. We demonstrate that solutions 
found previously are special cases of our general treatment. Finally in section 5 we 
briefly discuss the physical viability of our solutions. We emphasise that our simple 
approach of utilising the method of Frobenius for series yields a rich family of Einstein- 
Maxwell solutions in terms of elementary functions. This approach was utilised by 
Thirukkanesh and Maharaj— to produce electromagnetic solutions to the Einstein- 
Maxwell system that contains the Durgapal and Bannerji neutron star model—. 



2 The isotropic model 

On physical grounds it is necessary for the gravitational field to be static and spherically 
symmetric to describe the internal structure of a charged relativistic sphere. Therefore 
we assume that the interior of a spherically symmetric star is described by the line 
element 

ds^ = -e^-^^^de + e^^^'^dr^ + r\de^ + sin^ edcf?) (1) 

in Schwarzschild coordinates (t, r, 9, 0) where z/(r) and A(r) are arbitrary functions. For 
charged perfect fluids the Einstein-Maxwell system of field equations can be written in 
the form 

p + (2a) 

P - \e' (2b) 

P+\e^ (2c) 

y^e-\r^Ey (2d) 

for the line element ([T]). The field equations ([2]) are the same as those in Thirukkanesh 
and Maharaj—; we are utilising units in which the coupling constant and the speed of 
light are unity. The energy density p and the pressure p are measured relative to the 
comoving fluid 4-velocity = €""6^ and primes denote differentiation with respect to 
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the radial coordinate r. In the system ([2]), the quantities E and a are the electric field 
intensity and the proper charge density respectively. 

To integrate the system ([2]) it is necessary to choose two of the variables u, X, p, p 
or E. In our approach we specify A and E. In the integration procedure we make the 
choice 

where K is an arbitrary constant. The form ([3]) for the gravitational potential A re- 
stricts the geometry of the 3-dimensional hypersurfaces \t = constant} to be spheroidal. 
When K = the hypersurfaces {t = constant} become spherical. For the choice K = 
familiar spacetimes are regainable for particulars forms of the metric function e^^, e.g. 
the choice i/ = gives the metric of Einstein's universe. On eliminating p from ( ]2bl) 
and (l2cl) . for the particular form ([3]), we obtain 

\ r 

- {l-K){r/R')(^u'+^-^ + ^-^{l-KryR') (4) 

which is the condition of pressure isotropy with a nonzero electromagnetic field. 
It is convenient at this point to introduce the transformation 

^(a;)=e^M, = 1 - ^ (5) 
Then the condition of pressure isotropy (jlj) becomes 



:i - K + Kx^Y R^E^ 
x2 - 1 



[l~K + Kx^)i)-Kxi)+[- ^ + K{K = {] (6) 



in terms of new variables ip and x\ dots denote differentiation with respect to x. The 
Einstein-Maxwell system ([2]) implies 

l-K i3-K + Kx') 1 2 . 
P - i?2 (i-K + Kx^y'2 ^^^^ 

^ - R^il-LKx^) {-'i^^-')^l^' 

^ [2xE - {1 ~ x')E]^ 

R^{1 - x^){l ~ K + Kx^) ^ ' 

in terms of the variable x. Thus p, p, and o are defined in terms of in ([7]) 

The solution of the Einstein-Maxwell system depends on the integrability of ([6]). 
Clearly ([6]) is integrable once E is specified. A variety of choices for E is possible; 
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however only a few are physically reasonable. We need to choose E such that closed 
form solutions are possible. We make the choice 

2 ^ aKjx^ - 1) 

where a is constant. A similar form of E was also used by Sharma et aP and Tikekar 
and Singh—. The electric field intensity in ([8]) vanishes at the centre of the star, and 
remains continuous and bounded in the interior of the star for a wide range of values 
of the parameters a and K. Thus this choice for E is physically reasonable and useful 
in the study of the gravitational behaviour of charged stars. On substituting ([8]) into 
(E]) we obtain 

{1-K + Kx^)p^-Kx^ + K{a + K-I)^lj = (9) 

This is a second order differential equation which is linear in ip. We expect that our 
investigation of Q will produce viable models of charged stars since the special case 
a = yields models consistent with neutron stars. 



3 Series solution 

It is possible to express the solution of ^ in terms of special functions namely the 
Gegenbauer functions as demonstrated by Sharma et aP-. However that form of the 
solution is not particularly useful because of the analytic complexity of the special 
functions involved. In addition the role of parameters of physical interest, such as 
the spheroidal parameter K, is lost or obscured in the representation as Gegenbauer 
functions. The representation of the solutions in a simple form is necessary for a 
detailed physical analysis. Consequently we attempt to obtain a general solution to 
the differential equation ([9]) in a series form using the method of Frobenius. Later 
we will indicate that it is possible to extract solutions in terms of polynomials and 
algebraic functions for particular parameter values. 

As the point a; = is a regular point of ([2]), there exist two linearly independent 
solutions of the form of a power series with centre x = 0. Thus we assume 

oo 

^{x) = ^aix' (10) 

i=0 

where the constants are the coefficients of the series. For a legitimate solution we 
need to determine the coefficients Oj explicitly. On substituting (ITUi) into ^ we obtain 
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after simplification 

(1 - K)2.la2 + K{a + K- l)ao + [(1 - K)3.2a3 + K{a + K- 2)ai]x 

oo 

+ -K){i + + 2)ai+2 + K[a + K-1 + i{i - 2)]a,}x' = (11) 

i=2 

in increasing powers of x. For equation (fTTj) to be valid for all x in the interval of 
convergence we require 

{1-K)2.1a2 + K{a + K -l)ao = (12a) 

{1-K)3.2a3 + K{a + K -2)ai = (12b) 

{1- K){i + l){i + 2)ai+2 + K[a + K -l + i{i-2)]a, = 0,i>2 (12c) 

Equation (I12cl) is the the linear recurrence relation governing the structure of the 
solution. 

The recurrence relation (I12cl) consists of variable, rational coefficients. It does not 
fall in the known class of difference equations and has to be solved from first principles. 
It is possible to solve fll2cp using the principle of mathematical induction. We first 
consider the even coefficients Oq, 02, 04, ... Equation (I12al) imphes 

= (k^) + ^ ~ ^ + - 2)(2g - 4)]ao (13) 

where we have utilised the conventional symbol Y[ to denote multiplication for the first 
term. We now assume the pattern 



a2.= (^)'^n[« + ir-l + (2g-2)(2g-4)] 



ao (14) 



for the coefficient a2p which is the inductive step. We now establish that this is true 
for the next coefficient a2(p+i). Replacing i with 2p in f ll2cp we obtain 

K \ fa + K -l + 2p{2p-2] 
- yi^JK (2p+2)(2p+l) 1"'^ 
K \ fa + K -1 + 2p{2p - 2] 



X 



K-IJ \ (2p + 2)(2p + l 

I TTTTT I I [a + ^ - 1 + (2g - 2j(2g - 4jJao 



rT)'(^n[" + ^-l + ^2g-2)(2g-4)]c 

\ p+i p+i 

) [2{p + 1)]! + - ^ + - 2)(2g - 4)]ao (15) 



K-l 
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where we have used (1141) . Hence by mathematical induction all the even coefficients 
can be written in terms of the coefficient Qq. These coefficients generate a pattern 

""'^ = (k^) lit" + ^ - 1 + - 2)(2g - 4)]ao (16) 

for the even coefficients Oq, «45 • • • 

We can obtain a similar formula for the odd coefficients ai, 03, 05, . . . From fll2bp 
we have 

= (2 1 + 1)1 ri[« + ^ - 1 + (2g - l)(2g - 3)]ai (17) 

for the ffist term. We now assume 

a2p+i = ' (2p+i)! ri[« + ^ - 1 + (2g - l)(2g - 3)]ai (18) 

for the coefficient 02^+1. We then establish that this is true for the next coefficient 
a2(p+i)+i. Replacing i with (2p + 1) in (I12cp we obtain 

- ^T^ljy (2p + 3)(2p + 2) 

\ /a + ir- 1 + (2p+ l)(2p- 1) 



ir- 1; V (2p + 3)(2p + 2) 



(2^ + ^ - 1 + (2g - l)(2g - 3)]a, 

^ NP+I ^ P+1 



on utilising (ITS!) . Hence by using mathematical induction all the odd coefficients a2j+i 
can be written in terms of the coefficient ai. These coefficients generate a pattern which 
is clearly of the form 

a2.+i = (^-^ ^[a + K-l + (2g - l)(2g - 3)]ai (20) 

for the odd coefficients ai, 03, 05, . . . 

The coefficients a2i are generated from flTB]) . The coefficients a2j+i are generated 
from fl20l) . Hence the difference equation fll2cl) has been solved and all nonzero coef- 
ficients are expressible in terms of the leading coefficients oq and ai. From (fTOj) . (fT6l) 
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and (!20l) we establish that 



^{x) = ao|l + > J^r^) 7;^J|[a + K-l + (2g-2)(2g-4)]x 



i=l ^ ^ ^ ' g=l 



+ fli I a; + > , 1 ^7 — - ) 77T^^ I I [a + - 1 + (2g - l)(2g - ?>)\x^'+^ 

(21) 

where Oq and ai are arbitrary constants. Clearly (12T!) is of the form 

V'(a;) = ao^/'i(a;) + ai^2(a;) (22) 

where 

^i(^) = + (23a) 

(23b) 

are linearly independent solutions of (jH])- Thus we have found the general series solution 
to the differential equation ([9]) for the choice of the electromagnetic field E given in 
([8]). The solution (12T|) is expressed in terms of a series with real arguments unlike the 
complex arguments given by software packages. The series fl23ap and fl23bl) converge if 
there exist a nonnegative value for the radius of convergence. Note that the radius of 
convergence of the series is not less than the distance from the centre (x = 0) to the 
nearest root of the leading coefficient of the differential equation ([9]). Clearly this is 
possible for a wide range of values for K. 



4 Solutions with elementary functions 

It is interesting to observe that the series in (123!) terminates for restricted values of the 
parameters a and K. This will happen when a + K takes on specific integer values. 
Utilising this feature it is possible to generate solutions in terms of elementary functions 
by determining the specific restriction on a and K for a terminating series. Solutions 
in terms of polynomials and algebraic functions can be found. We use the recurrence 
relation (I12cp . rather than the series fl23|) . to find the elementary solutions as this is 
simpler. 
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4.1 Polynomial solutions 

We first consider polynomials of even degree. It is convenient to set 

I = 2(j-l) (24a) 
K + a = 2-{2n-l)'^ (24b) 

where n > 1 is fixed integer in f ll2cl) . This leads to 

{n + j-2){n-j + l) 
= -7 2j{2j-l) "^^-^ ^^^^ 

where we have set 7 = 4 — ^n{n-i)+a - ^ote that (12^ implies a2(n+i) = 0. Conse- 
quently the remaining coefficients a2(n+2)5 (i2{n+3), ■ ■ ■ vanish. Equation fl2^ may be 
solved to yield 

= (-')'t^7)W' '''' 

where we have set Oq = ^(n-i) ■ With the help of (l26ll we can express the polynomial 
in even powers of x in the form 

^ (n-j)!(2j)! 

for is: + a = 2 - (2n - 1)^. 

We now consider polynomials of odd degree. For this case we let 

i = 2(j-l) + l (28a) 
K + a = 2(1 -2n^) (28b) 

where > is fixed integer in fll2cl) . We obtain 

{n+j-l){n-j + l) 
a2,+i = -/X 2j(2j + l) ^^^^ 

where we have set /x = 4 — -(J^^t^y^- We observe that (129|) implies a2(n+i)+i = 0. 
Consequently the remaining coefficients a2{n+2)+i? a'2(n+3)+i5 • • • vanish. Equation 
(12^ can be solved to yield 

(n A- j — l)\ 

°^^^' = '-^)' (n-,)K2, + l)! - (3°) 

where we have set ai = ^. With the assistance of ( l30l) we can express the polynomial 
in odd powers of x as 

«.W-g(-.)'(;r^^Jyg^^-- (31) 

for K + a = 2(1 -2n2). 

The polynomial solutions (1271) and (I5T1) comprise the first solution of for appro- 
priate values of -|- a. 
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4.2 Algebraic solutions 

We take the second solution of (jH]) to be of the form 

ij{x) = u{x){l - K + Kx^f'^ (32) 

when u{x) is an arbitrary polynomiaL Particular solutions found in the past are special 
cases of this general form; the factor (1 — K + Kx'^Y^'^ helps to simplify the integration 
process. This motivates the algebraic form for ip as a. generic solution to the differential 
equation ([9]). On substituting in ([9]) we obtain after simplification 

{l-K + Kx^)j^ + 5Kx^ + K{a + K + 2)u = (33) 

which is linear differential equation for u{x). 

As in §4.1 we can find two classes of polynomial solutions for u{x), in even powers 
of X and in odd powers of x, for certain values oi K + a. As the point x = is a regular 
point of (l33l) . there exists two linearly independent solutions of the form of power series 
with centre x = 0. Therefore we can write 



u[x} = ^bix' (34) 

i=0 

where 6j are the coefficients of the series. Substituting flM|) in fl55]) we obtain 
(1 - K)2.1b2 + K{a + K + 2)bo + [(1 - K)3.2b3 + K{a + K + 7)bi]x 

oo 

+ ^{(1 -K){i + 2){i + l)k+2 + K[a + K + 2 + i{i + A)]bi}x' = (35) 

i=2 

For equation (!35|) to hold true for all x we require that 

(1- K)2.1b2 + K{a + K + 2)bo = (36a) 
{1- K)3.2bs + K{a + K + 7)bi = (36b) 
{l-K){t + 2){t + l)bi+2 + K[a + K + 2 + i{i + 4)]bi = 0,i>2 (36c) 

which governs the coefficients. 

We first consider even powers of x. Replacing i with 2(j — 1) and assuming K + a = 
2(1 — 2n^) in fl36cp . where n > is fixed integer, we obtain 

{n + j){n-j) 
= 2j(2j-l) ^^^^ 

where we have set = 4 — ^^2^i^^ - From ( 1371) we have that 62™ = and subsequent 
coefficients &2(n+i)5 &2(n+2)5 • ■ • vanish. Then (1371) has the solution 



(n-j-l)!(2j)!' 



b2j = (-/i) \„ \. < J < n - 1 (3^ 
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where we have set 60 = On using (IMIl and (!38l) the polynomial in even powers of x 
leads to the expression 

= (1 - + ^•-^)"^ i:'-^) V- j-l)!(2,)! "'' 

for + a = 2(1 - 2n^). 

We now consider odd powers of x. Replacing i with 2(j — 1) + 1 and assuming 
K + a = 2 — {2n — 1)^ in fl36cp . where n > 1 is fixed integer, we obtain 

= -7 2j(2j + l) 

where we have set 7 = 4 — ^n{n-i)+a - ^^^m (HOl) we have that 62{n-i)+i = and 
subsequent coefficients 62n+i, &2(n+i)+i, • • • vanish. Then equation pUj) has the solution 



2.+1 = (-7) \,. / J/;, , < J < n - 2 (41) 



(n + j)! 
(n-j-2)!(2j + l)! 

where we have set bi = n{n — 1). On using fl34l) and fHT]) the polynomial in odd powers 
of X leads to the result 



n-2 



Ux) = {1-K + Kxy'^Y.^-^yi t^uo■^^^^ ''"'^" ^^2) 

^ (n-j -2)!(2j + l)! 

for + a = 2 - (2n - 1)^. 

The algebraic solutions fl5^ and fH21) comprise the second solution of (jH]) for ap- 
propriate values of + a. The solutions fl39|) and fj^2l) are expressed as products of 
algebraic functions and polynomials, and they are clearly linearly independent from 
(13T1) and fl271) . respectively. 



4.3 Elementary functions 

We have obtained two classes of polynomial solutions (1271) and (!3T|) in §4.1 to the 
differential equation Also we have found two classes of algebraic solutions fl5^ 
and (j^2|) in §4.2. By collecting these results we can express the general solution to ([9]) 
in two categories. The first category of solution for il){x) = f{x) is given by 

fix) = Af,{x) + BUx) 



{n-j)\{2jy. 



x^^ 



n— 2 / .s^i 

+ + (43) 
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for the values 

4 

7 = 4- — (44a) 

4n{n - 1) + a ^ ^ 

K + a = 2-(2n-l)2. (44b) 
The second category of solution for iIj{x) = g{x) has the form 
g{x) = Agi{x) + Bg2{x) 



n 



[n-j)\{2j + l)\ 
B(l-K + K.-f" U-^y ^^:';^,^,^/ ' (45) 



for the values 



= 4- . 2 ^ (46a) 
ir + a = 2(1 -2n2) (46b) 

where A and -B are arbitrary constants. 

It is remarkable that these solutions are expressed completely as combinations of 
polynomials and algebraic functions. From our general class of solutions (H3|l and (l45l) 
it is possible to generate particular solutions found previously. Consider one example 
with a = and K = — 7(n = 2). Then 7 = | and it is easy to verify that equation (l43l) 
becomes ^ 

ij = A(^l-L' + ^x'^ +Bx(^l- Ix'^ ' . (47) 

Thus we have regained the Tikekar- solution for a superdense neutron star from our 
general solutions. Many other particular solutions found in the literature are also 
contained in our general solutions, eg. the model of Patel and Koppar— . The solutions 
(H3l) and (H5l) reduce to the Maharaj and Leach^ model when a = 0. Our solutions 
are applicable to a charged superdense star with spheroidal geometry. When a = we 
obtain uncharged relativistic stars which model ultradense barotropic matter. 



5 Physical conditions 

In the general solution of Sharma et al^ it is not possible to isolate the spheroidal 
parameter K as that solution is given in terms of special functions. Our solutions are 
in terms of simple elementary functions which facilitate a study of the physical features, 
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in particular the role of K. The exact solutions (H3!) and (H5!) make it possible to analyse 
the role of the spheroidal parameter K and it's connection to the electromagnetic field. 
In particular it is possible to make the following comment about the special role that 
the spheroidal parameter K has in charged solutions. The form of the solution for the 
uncharged relativistic star is similar to fH3l) and fHSl) : however the models are different 
because the coefficients of the polynomials (namely 7 and /x) differ by the parameter 
a. If the parameter a > then we observe that 

-ft'(a^O) < -^(a=0) 

Hence the presence of charge directly affects the spheroidal geometry through the 
parameter K. The geometry of the hypersurfaces {t = constant} in the spacetime 
manifold is related to the electromagnetic field via the relationships K = 2— (2n— 1)^— a 
and K = 2(1 — 2n^) — a. Such explicit relationships connecting the spacetime geometry 
to the energy momentum (or electromagnetic field) are rare in exact solutions. The 
presence of charge a decreases the value of the spheroidal parameter K in our solutions. 

We make a few brief comments about the physics of the models found in this paper. 
If < < l(a < 0) then p remains positive in the region 

which restricts the size of the configuration. When K < 0(a > 0) there is no restriction 
on p. Hence p is positive in the interior of the star. It is clear from fl7a|) and ([8]) that 
dp/dr < ioT K < 0(a > 0). Consequently the energy density decreases from the 
centre to the boundary. For the pressure to be vanish at the boundary r = a we 
require that 



1 \ 1 - aVi?2 1 aKa^/R^ , , 

' \ \ = f49) 



where t/' is given by ( |43|) or ( l45l) . This will constrain the values of the constants A 
and B. The solution of the Einstein-Maxwell equations for r > a is given by the 
Reissner-Nordstrom metric as 

= - (1 - ^ + ^) + (1 - ^ + ^) ' + + sin' ^#'), (50) 

where m and q are the total mass and charge of the star respectively. To match the 
line element ([T]) with the Reissner-Nordstrom metric ( |50l) across the boundary at r = a 
we require the continuity of the gravitational potentials and of the radial electric field 
at r = a. Continuity of the gravitational potentials yields the relationships between 
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the constants A, B, K, a and R as 



a a? 



[Ai)i{a) + Bi)2{a)Y (51a) 



2m q\a)\-' 1 - Ka'/R' 
'"V + ^J = 1 - aVR- ^''^^ 

The continuity of electric field yields the form 

for the charge at the boundary. This shows that continuity of the metric functions 
across the boundary r = a is easily achieved. The matching conditions at r = a may 
place restrictions on the metric coefficients v and its first derivative for uncharged mat- 
ter; and the pressure may be nonzero if there is a surface layer of charge. However 
there are sufficient free parameters to satisfy the necessary conditions that arises for a 
particular spheroidal model. It is interesting to note that our solutions may be inter- 
preted as models for relativistic anisotropic stars where the parameter a plays a role 
of the anisotropy factor. Isotropic and uncharged stars can be regained when a = 0. 
Chaisi and Maharaj— , Dev and Gleiser— and Maharaj and Chaisi— provide some 
recent treatments involving the physics of anisotropic matter. 
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